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REAL-TIME SIMULATION OF AN AUTOMOTIVE GAS TURBINE USING 


THE HYBRID COMPUTER 

William Costakls and Walter C. Merrill 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


SUMMARY 

A hybrid computer simulation of an Advanced Automotive Gas Turbine Power- 
train System Is reported. The system consists of a gas turbine engine, an 
automotive drlvetraln with four-speed automatic transmission, and a control 
system. Generally, dynamic performance Is simulated on the analog portion of 
the hybrid computer while most of the steady-state performance characteristics 
are calculated to run faster than real-time and makes this simulation a useful 
tool for a variety of analytical studies. 


INTRODUCTION 

This report describes a hybrid computer simulation of an Advanced Automo- 
tive Gas Turbine (AGT) Powertrain System, designated the AGT101. The AGT101 
Is being developed by the Garrett Turbine Engine Company under NASA contract 
DEN3-167 sponsored by the Department of Energy (ref. 1). The AGT101 repre- 
sents a significant advancement In gas turbine state-of-the-art. The engine 
Incorporates ceramic parts to enhance high temperature operation and to pro- 
vide low cost manufacturing potential for automotive application. 

As part of the development of the AGT101 a digital computer simulation 
was developed to Investigate various control strategies and to study engine/ 
drlvetraln control dynamic performance. However, since a real-time simulation 
capability has proven to be very useful and productive tool for controls 
research, a hybrid computer simulation Implementation of the AGT101 was 
desired. Therefore the digital simulation was converted to a hybrid computer 
based simulation and Is reported herein. 

The model, as defined by the engine manufacturer's digital simulation, 
consists of three basic subsystems: (1) the engine, (2) the drlvetraln, and 

(3) the control. Since the hybrid simulation Is to serve as a controls 
research testbed, the detailed control subsystem supplied by the manufacturer 
was not required. However, to facilitate the comparison between the hybrid 

and digital simulations, the control subsystem was Included In the hybrid 
simulation but In a simplified form. These simplifications result In some 
minor discrepancies between the hybrid and digital simulations, particularly 
In small transients close to Idle speed. However, these discrepancies do not 
seriously Impact the usefulness of the hybrid simulation. Further, dynamic 
and steady-state comparisons between the hybrid and digital simulations of the 
engine and drlvetraln components are good. These comparisons Indicate that 
the hybrid simulation can be used as a realistic controls testbed facility. 


The remainder of the report describes the engine, drlvetraln, and control 
model, the actual hybrid Implementation and the results obtained from the 
hybrid simulation. 


SYMBOLS 

ALPO torque converter output angular acceleration 

ALPOUT gearbox output angular acceleration 

ALPS1 engine angular acceleration 

BETA VSTC vane angle, deg 

CD aerodynamic drag coefficient 

CSLIP clutch slip, rpm 

ONE speed error, ANE/NE 

DNERATE rate of change of speed set point 

DTQ/DT4 engine torque correction factor 

DT5.1D4 hot side of regenerator temperature correction factor 

DWDT4 air flow rate correction factor 

EFAXLE rear axle gear efficiency 

EFREG regenerator effectiveness 

ENGS normalized engine speed 

EPS.ETR gear box efficiency 

FACC force available for acceleration 

GAMMA throttle position, deg 

GAMI throttle position Integral control, deg 

GAMP throttle position proportional control, deg 

GRAXL rear axle gear ratio 

GRTR transmission gear ratio 

GRZ gear ratio of gear between gearbox and transmission 

HPAC power consumption by accessories, hp 

HPACC power available for acceleration, hp 

HPCAR power at Input to gearbox, hp 

HPDS power of the drlveshaft, hp 

HPLOSS transmission power loss, hp 

HPOTR power at output of gearbox, hp 

HPRL power required to overcome road load, hp 

HPWHL available power at the wheels, hp 

HV fuel heating valve 

IAC Inertia of accessories, lb-ft-sec 2 

IGV compressor Inlet guide vane position 

I G VI component of IGV which Is a function of engine speed error 

IGV2 component of IGV which Is a function of GAM 

I PI 2 Inertia of planet gear 1 and 2 

IP3 Inertia of planet gear 3 

IR1 Inertia transformed to ring gear 1 

IR1 ' Inertia of ring gear 1 

IR2 Inertia of ring gear 2 

IS1 Inertia of sun gear 1 

ITQC Inertia of Input side of torque converter 

ITQU Inertia of output side of torque converter and Input side of clutch 

ITQO Inertia of outside clutch and planet gears 1 and 2 

JE Inertia of engine up to torque converter, lb-ft-sec 2 

KPB4 VSTC vane control gain 

KSET engine speed set point gain 
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NC 

NCTR 

NOS 

NE 

NECOM 

NEP 

NESET 

NOTR 

NP1 

NP2 

NP3 

NR1 

NR2 

NS1 

NS2 

NWHL 

REVMI 

ROL1 

RN 

SIGMA(o) 

SPIN1 

SPIN2 

SPRTIO 

TCKI 

TCSR 

TCTR 

TQAC 

TQACC 

TQCAR 

TQCTR 

TQDS 

TQE 

TQOTR 

TQRL 

TQT 

TQTOT 

TQWHL 

T3.1 

T3.6 

T4 

T4ERR 

T4SET 

T5.1 

T5.1ERR 

T5.1SET 

VEL 

VELC 

VSTC 

WAIR 

WAWF 

WF 

WO 


combustion efficiency 

speed of torque converter Input, rpm 

drlveshaft speed, rpm 

engine speed, rpm 

engine speed command, percent 

engine speed, percent 

engine speed set point, percent 

gearbox output speed, rpm 

number of teeth In planet gear 1 

number of teeth In planet gear 2 

number of teeth In planet gear 3 

number of teeth In ring gear 1 

number of teeth In ring gear 2 

number of teeth In sun gear 1 

number of teeth In sun gear 2 

wheel speed, rpm 

revolution per mile of tire 

rolling friction factor 

normal vehicle rolling resistance 

drag coefficient gain on factor 

ambient transmission power loss 

transmission power loss factor 

ration of torque converter Input speed to accessory design speed 

torque converter Input capacity factor, rpm/ /FI 

torque converter speed ratio 

torque converter temperature ratio 

accessory torque, ft-lb 

torque available to accelerate vehicle, ft-lb 

torque Input to gearbox, ft-lb 

torque at Input of torque converter 

torque on drlveshaft, ft-lb 

unbalance torque accelerating engine, ft-lb 

gearbox output torque, ft-lb 

torque required to overcome road load, ft-lb 

uncorrected engine torque 

engine torque, ft-lb 

torque to wheels, ft-lb 

regenerator coldslde Inlet temperature, °R 

combustor Inlet temperature, °R 

turbine Inlet temperature, °R 

difference between turbine Inlet temperature and Its set temperature 
point (T4SET), °R 
T4 set temperature point, 2960° R 
regenerator hot side Inlet temperature, °R 

difference between regenerator hot side temperature and Its set 
temperature point, ®R 
T5.1 set temperature point, 2460° R 
vehicle speed 
vehicle speed command 
variable stator torque converter 
airflow, lb/sec 

combustor gas flow rate, lb/sec 
fuel flow rate, lb/hr 

Inertia of transmission through wheels, ft-lb 
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WO speed of torque converter output, rpm 

WT vehicle weight, lb 

WTEQ equivalent vehicle weight, lb 

WFGAIN governor gain, lb/hr/NEP 


SYSTEM MODEL DESCRIPTION 

This section describes the model of AGT101 system. The model Includes a 
description of the engine, drlvetraln, and control subsystem. 


Engine 

The engine, shown In figure 1, Is a regenerated single-shaft gas turbine 
engine which Is flat-rated at 100 horsepower. The engine consists of a 
single-stage, 5 to 1 pressure ratio, backward-swept centrifugal compressor and 
single-stage ceramic radial Inflow turbine mounted on a common shaft. Maximum 
rotor speed for this compressor and turbine configuration Is 100 000 rpm. 
Compressor Inlet guide vanes are used to vary engine airflow, and combustor 
fuel flow Is used to vary combustor temperature. 

The engine Is modeled as shown In figure 2. Model equations are pre- 
sented In the appendix. The engine model Is composed of three parts, the 
engine performance maps, the regenerator, and the combustor. Engine airflow, 
torque, and regenerator hot-side and cold-side Inlet temperatures are deter- 
mined as dependent variables from engine performance data maps using turbine 
Inlet temperature, rotor speed, and Inlet guide vane positions as Independent 
variables. Combustor Inlet (regenerator exhaust) temperature Is determined by 
an energy balance In the regenerator which Incorporates the efficiency of heat 
transfer In the regenerator. Regenerator heat storage Is modeled with a 
regenerator thermal time constant which Is dependent upon engine airflow. 
Combustor operation Is modeled by 'curve fit 1 approximations to the perfor- 
mance of an analytically based, detailed thermodynamic combustor simulation. 
This simplified combustor model generates the combustor exhaust (turbine 
Inlet) temperature from combustor Inlet temperature, air flow, and fuel flow. 
Finally, engine speed Is determined by a conservation of angular momentum 
equation. The difference between total developed engine torque and the torque 
required by the drlvetraln Is used to accelerate the engine rotor. 


Drlvetraln 

For the purposes of this report the drlvetraln subsystem Includes the 
gearbox, the transmission, and Important physical characteristics of the 
vehicle. The gearbox Is a split path differential design accepting power 
directly from the engine. Power Is split In the differential planetary gears 
with a portion going directly to the transmission. The remaining power goes 
to the variable stator torque converter (VSTC). The VSTC provides a variable 
speed ratio output that Is fed back to the planetary carrier and combines with 
the power delivered to the transmission. The transmission Is the Ford four- 
speed Automatic Overdrive (AOD) production transmission. The transmission 
model consists of four forward gears. Reverse gearing Is not considered In 
this simulation. Shifting logic Is a function of a vehicle velocity. The 
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transmission couples the output power of the gearbox with the vehicle. The 
vehicle drlvetraln subsystem Is composed of the drlveshaft, the differential, 
the rear axle of the car, and the wheels. Vehicle weight, Inertia, and road 
load characteristics are also Included In the simulation. The subject vehicle 
of this simulation Is a 3000-pound vehicle of the Ford Fairmont class. 

The drlvetraln model Is also shown In figure 2 with model equations given 
In the appendix. Basic Inputs to the drlvetraln are the VSTC control vari- 
able, BETA, and engine rotor speed. The outputs of the drlvetraln model are 
vehicle required torque as well as Important Internal drlvetraln variables 
such as vehicle velocity. 

The VSTC vane angle, BETA, along with a calculated VSTC speed ratio are 
used as Inputs to performance maps (obtained foom test data) of the VSTC. The 
outputs are VSTC torque ratio and VSTC capacity factor. These outputs along 
with efficiency of the gearbox, transmission, and axle, the gear ratio of the 
transmission, vehicle velocity, and accessory load are used to calculate the 
vehicle required torque. Vehicle velocity Is calculated from a conservation 
of angular momentum equation. The difference between the torque applied to 
the wheel and the torque due to the road load accelerates the vehicle. Road 
load Is found as a function of vehicle weight, velocity, and grade of the road. 


Control 

The control model Is shown schematically In figure 2 and the model equa- 
tions are given In the appendix. Basically, the control has three Independent 
or manipulated variables, compressor Inlet guide vane position (IGV), fuel 
flow rate (WF), and VSTC vane angle (BETA). Inlet guide vane position Is used 
to optimize compressor performance with respect to shaft speed and Is sched- 
uled open-loop as a function of commanded rotor speed error and the throttle 
pedal angle (GAMMA). The fuel flow rate Is determined by a closed-loop 
proportional plus Integral control on engine rotor speed. Command engine 
speed Is a direct function of GAMMA. Finally, VSTC vane angle (BETA) Is 
determined by closed-loop proportional control on either combustor exhaust 
temperature or turbine exhaust temperature. One, or the other, of these two 
temperatures Is controlled to Its respective constant, predetermined setpoint 
value. The maximum of the two temperature errors (setpoint minus actual) 
determines which loop Is selected as the active control loop. 


HYBRID SIMULATION 

The equations describing the AGT101 power train system were Implemented 
on the Lewis Research Center's PACER600 hybrid computing system. This system 
consists of an EAI PACER 100 digital computer, a model 680 analog computer and 
a model 681 analog computer. 

The digital computer was used primarily to perform the generation of the 
engine's performance maps and the torque converter maps. The digital computer 
was also used to Implement the shifting logic and to calculate the ambient 
transmission power loss (SPIN1) and the transmission power loss factor 
(SPIN2). The scaled-fraction variables and the digital function generator 
programs MAP and MAPL, used In the digital programs, are discussed In 
reference 2. 
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The remaining calculations were performed on the analog computer. The 
analog computer performed such computations as summing, Integration with 
respect to time, multiplication, unlvarlable function generation, and so 
forth, which are suitable to analog computing. Strip-chart recorders were 
used to continuously monitor and record the computed variables. 


RESULTS AND DISCUSSION 

The Intent of the hybrid simulation was to demonstrate the feasibility of 
using the hybrid computer to simulate. In real-time, AGT101 performance. The 
simulation then could be used as a test bed for evaluating various research 
control strategies. 

The basis for the hybrid simulation was the digital simulation provided 
by the engine manufacturer. Some digital simulation results are presented for 
comparison with hybrid simulation results to establish validity of the hybrid 
results . 

A comparison between steady-state hybrid and digital simulation results 
Is presented In Table I. Results for throttle settings of 20°, 40° and 50° 
are shown In Table I. Steady-state results for both the digital and hybrid 
simulations were obtained by first Increasing the throttle from 0° to the 
final desired throttle setting In a step-wise fashion, and then allowing the 
simulation sufficient time to reach Its steady-state point. 

Twenty-three engine variables are compared In Table I. For every one of 
these variables the hybrid results compared well with the digital results. 

The transient data are responses of selected engine variables to a step 
throttle command under zero road grade conditions. Prior to the throttle step 
Input, the simulation was set at the Idle (GAMMA = 0) Initial conditions. 
Transient data from the hybrid simulations were obtained for throttle command 
step Inputs from 0° to 20° and 40°, with the results presented In figures 3 
and 4. Transient responses of eleven engine variables, using the digital 
simulation, are presented In figures 5 and 6. The responses are for throttle 
step Inputs of 20° and 40°. A comparison between the hybrid and the digital 
simulation dynamic responses shows that they are quite similar. Figures 3 and 
5 which are the results of the 20° throttle step Input for the hybrid and 
digital simulations, respectively. Indicate that the transmission shifted 
through the first three gears only, and the vehicle attained Its final steady- 
state velocity In the third gear. The results for the 40° throttle step 
Input, shown In figure 4 for the hybrid simulation, show that the transmission 
shifts through all four gears and the vehicle attains Its final steady-state 
velocity In the fourth gear. The digital simulation results shown In figure 6 
are similar. 

Also the results In figures 3 and 4 Indicate that the smaller the throt- 
tle step Input the longer It takes to shift from one gear to the next. Con- 
sequently, the larger the step Input the faster the vehicle attained Its final 
steady-state velocity. The digital simulation results showed the same trend. 

The only difference of any significance between the hybrid and the 
digital simulation with regards to the transmission was In the shift point 
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timing. The hybrid simulation used somewhat less time between shift points 
than the digital simulation. This difference, however, was expected since the 
actuator dynamics associated with the IGV's and VSTC's were not Included In 
the hybrid simulation. 

The other area where the simulation differed was at an Idle or near Idle 
condition. In the digital simulation with the transmission In first gear, 
throttle command at zero degrees and zero road grade conditions, the vehicle 
attained a steady speed of 0.26 miles per hour. Prior to the throttle command 
step Input In the hybrid simulation, the vehicle attained a steady velocity of 
15.68 miles per hour. This can be explained by the fact that hybrid simula- 
tion did not Include the detailed control employed In the digital simulation. 
As an example, the hybrid simulation did not provide for the reduction of fuel 
flow when T4 and T5.1 exceeded the over temperature limits. It also did not 
provide for Integration of IGV's to open position If either T5.1 or T4 
exceeded their set points. The net result was that the hybrid simulation 
would run richer at Idle and consequently would stabilize at a higher vehicle 
velocity than the digital. Since this was a control system related discrep- 
ancy, It has no significant effect on the evaluation of the hybrid simulation. 

The dynamic response of all the torques, the position of the VSTC vanes 
(BETA), percent engine speed (NEP), vehicle velocity (VEL) and fuel flow (WFL) 
to a 20° and 40° throttle command shown In figures 3 and 4, respectively, com- 
pared well with the digital simulation results shown In figures 5 and 6. 

The hybrid simulation compressor Inlet guide vane (IGV) position dynamic 
response shown In figures 3 and 4 varied somewhat from the digital simulation 
response shown In figures 5(b) and 6(b). However, since there was significant 
control logic In the digital simulation that affected the IGV's directly which 
did not exist In the hybrid, this small difference In the IGV response was not 
considered Important. 


SUMMARY OF RESULTS 

This report describes a real-time hybrid computer simulation of the 
AGT101 automatic gas turbine powertrain system. The AGT101 real-time simula- 
tion was modeled after the manufacturer's digital simulation. Since the real- 
time simulation was Intended to serve as a control research testbed, the 
detailed control supplied by the manufacturer's digital simulation was not 
Included. Steady-state data are presented for three throttle positions at 
zero road grade conditions. Transient data are presented In the form of 
engine variable responses to two throttle step Inputs under zero road grade 
conditions. Both steady state and transient data of the real-time simulation 
compared well with the digital simulation data. Differences between the two 
sets of data at Idle or near Idle condition were due to the exclusion of the 
control supplied with the digital simulations from the real-time simulation. 
The real-time simulation shifted somewhat faster than the digital simulation 
because the actuator dynamics associated with the IGV's and VSTC vanes were 
not Included In the simulation. 

The simulation structure described In this report resulted In stable 
real-time operation. 
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APPENDIX - EQUATIONS 


I. Simulated Driver Commands: 

GAMMA = GAM I + GAMP 
GAMP = KV1 AVEL 
GAMI = KV2 AVEL dt 
AVEL = VELC - VEL 

II. IGV Control: 

IGV = IGV1 + IGV2 

IGV1 = f(DNE) = -ONE 

IGV2 = f( GAMMA) 

Tru . IGV 

' 1 * t, gv S IGV 

III. Fuel Speed Control: 

NECOM = f( GAMMA) 

WF = WFGAIN(NESET - NEP) 
KSET 


NESET = 


(DNERATE) 


DNERATE = f(DNE) 
NECOM - NEP 


DNE = 
WFL = 


NEP 

WF 


1 + T p S 


O.Ol 


IV. Engine Maps 

Engine air flow: 

WAIRM = f ( NEP, IGVL) 

WAIRC = WAIRM - DWDT4(T4 
7.5 IGVL 


DWDT4 = 


WAIRL = 


10 3 10 ' 

WAIRC 

1 + t.S 
A 


t A = 


= 0.01 

WFGAIN = 3 
KSET = 4 


- 2960) 

1 .01 
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Engine torque: 


TQT = f (NEP, IGVL) 
DTQDT4 = f (NEP, IGVL) 


TQTOT = TQT + DTQDT4 


T. - 2760 
4 

1000 


- 5xlO~ 5 


(T 4 - 2760) 


Regenerator temperatures: 


T3.1 = f (NEP, IGVL) 

T5.1H = f (NEP, IGVL) 

DT5.1D4 = 0.77 + 0.006170NEP - ( 0 . 61 3642x1 0“ 4 NEP 2 ) 
T5.1 = T5.1M + DT5.1D4(T4 - 2960) 


V. 


Regenerator 

EFREG = 0.98525 - 0.01074WAIRL - 0.1004WAIRL 2 
T3.6 = T3.1 + EFREG(T5.1 - T3.1) 

T3.6 


T3.6L 


t REG 


1 * t reg s 

4.25 


= 0.33 + 


WAIRL 


0.03617WAIRL 3 


VI. Combustor 


T4 _ r(ECOMB) 2 + DC0MB1 0, 5 - ECOMB 
0. 224x1 0 4 WAWF 


WAWF = WAIRL + 


ACOMB = 


1 


0.664 


WFL 

3600 

10,65 + 3600WAIRL^ 


BCOMB = WAIRL (0.24 + 0.61xl0~ 6 Y) (T3.6L - 536.7) 
CCOMB = n c (HV) ^ 


DCOMB = 0.448x1 O' 4 WAWF 
ECOMB = 0.2355 WAWF 


BCOMB + CCOMB + 129.3WAWF 
ACOMB 
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Y = (T3.6L ) 1 - 35 
VII. VSTC Vane Control 


T5.1L = 


T5.1 

1 + t th S 


t TH 


= 4 


1 + T.S T - ) 

T5.1LL = 7 — — c T5.1L 
• + T oS 

2 t2 

T5.1ERR = T5.1LL - T5.1SET 
T4 


T4L = 


1 * , TH S 


= 4 

= 0.8 

T5.1SET = 2460° 


R 


1 + t S 

T4LL = 7 — — ^ T4L 
I + 

T4ERR = T4LL - T4SET T4SET = 2960° R 
BETA = 84 + KPB4(TERR) KPB4 = 0.8 
TERR = Max(T5.1 ERR.T4ERR) 


VIII. Gear Box System 

NOTR = NDS( GRTR) (GRZ) GRZ = 1.0405 
NDS = NWHL(GRAXL) GRAXL = 2.730 
REVMI 


NWHL = 


60 

GRTR = f (GEAR) 


(VEL) 


HPOTR = 


TQOTR(NOTR) 


5252 


TQOTR = E 6 


( F 6 ) (TQCAR) - G 6 ( ALPS1 ) + IP12(ALP0UT) 

+ TQCTR(TCTR) 


m] 


EPS = ETR 


ALPS1 


2ir dNE 
60 dt 


2ir(12)(9.5493) 
60( JE) 


TQE 


ALPOUT = ^ 


REVMI 

60 


( GRAXL) ( GRTR) ( GRZ) 


dVEL 

dt 


ALPO(ITQU) 
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A6 


■(») 

[ NS1 ( NP2) 1 
NR1 (NP1 ) J 


(IP12)(EPS) 2 = — 3 4~~ EPS 2 


10 ' 


pc NS1(NP2) n n , flifl 

66 = 1 NRl(NPl) 1 ' °* 06948 

NPJ.(NRl) 

C6 1 NS1 (NP2) 15,8571 


D6 _ f NS2(NS1)(NP2) 


NR1 (NP1 )(NR2) 



0.05256 

(EPS) 6 


E6 _ NPl ( NS1 ) _ q 12771 

tb ‘ NS1(NP2) " u - uin 

F6 = NP2( EPS) + (NP1)(EPS) 2 = (42 + 96EPS)EPS 


G6 = NPl ( IS1 )( EPS) * NP2( IS1 ) EPS + 


/ 1 PI 2 \ 

Inpi ) 


(NS1) = 0 


where 

IS1 = IP3 = IAC = IR2 = 0 
IP1 2 = 0.0004416 
IR1 ' = 0.036583 

2 


ITQC = 0.075 
ITQU = 0.0104166 
ITQO = 0.008125 


IR1 = IR1 ' + IP3 
and 


(hsz\ 

\ND3 J 


+ (IR2 + IAC + ITQC) 


( NS2 \ 
y ND2 J 


NS1 = 22 NPl = 96 NP2 = 42 NR1 = 143 

NS2 = 171 NP3 = 23 NR2 = 219 


Clutch Slip: 

For NEP > 56 or GAH > 1 

CSLIP . [a|r(NE) - ((As lf f i ii)) (NOTR)] TCSR - 

NPl ( NPl )(NR2) 

M;>K = NS1 (NP2)(NS2) 

BSR = NS1 ( NP2) 

CSR = (NS1 )(NP2) + (NPl ) ( NR1 ) 


.00229 


2 


NOTR 
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For NEP < 56 or GAM < 1 
CSLIP = 0 


Efficiency: 

ETR = AGE + BGE(HPCAR) + CGE(HPCAR) 2 + DGE(HPCAR) 3 
HPCAR - ISffSp 

AGE = 0.97015 + ENGS(0. 0307903 - 0.07365 ENGS) 

BGE = 5.2315xl0“4 + ENGS( 2. 1 5x1 0~ 4 + 1.0562xl0" 3 ENGS) 

CGE = 3.68359xl0“k _ ENGS(1 .20716 x 1 0“ 5 + 2.14907xl0- 6 ENGS) 
DGE = 7.47466xl0" 9 + ENGS(6.856xlQ- 8 - 1 . 5756x1 0~ 8 ENGS) 
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Torque Converter: 

ASR(BSR)(NOTR + CSLIPL) 


TCSR = 


NCTR = 


BSR(NE) - NOTR(CSR) 
WO 


TCSR 

WO = NOTR + CSLIPL 
TCKI = f (TCSR, BETA) 
TQCTR = ( NCTR/TCKI) 2 
TCTR = f (TCSR, BETA) 

HPCTR = KMggn 


IX. Accessory Torque and Power Consumption 
TQAC = 5252( HPAC/NCTR) 

HPAC = 2 .342( SPRTIO) 1 - 01 38 
SPRTIO = MAX(0. 001 .SPRTIO') 

SPRTIO' = NCTR/4000 


X. Engine Acceleration Torque and Speed 
TQE = TQTOT - TQCAR 
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dNE 114.6 
dt = JE 

NE 


TQE JE = 0.0052 


NEP = 


10 ^ 


= jlSl + A6 + 86 ALPS1 

- [ A6 * B6(C6) ffl " 


ALPOUT + ^TQAC + TQCTR^ 06 


) 


XI. Transmission Efficiency 


HPLOSS = SPIN! + SPIN2(TQDS) 
HPLOSS = Hax(0, HPLOSS) 


For gear 1 : 

SPIN1 = 0.09459 - 0. 08389xl0~ 3 NDS + 5 . 1 68x1 0“ 6 NDS 2 

- 4.020x10-9 nds 3 + 1.295X10" 12 NDS 4 

SPIN2 = 0.0003162 + 0.0247x10" 3 NDS - 0.02772xl0" 6 NDS 2 

+ 0.01043X10- 9 NDS 3 - 0.0002181X10" 12 NDS 4 


For gear 2: 

SPIN1 = 0.05734 + 0.4648x10-3 NDS + 0.7247xl0" 6 NDS 2 

- 0.1023x10-9 NDS 3 + 0.01 740xl0- 12 NDS 4 

SPIN2 = 0.0004343 + 0-01 31 Oxl 0" 3 NDS - 0.009127X10" 6 NDS 2 

+ 0.002259x10-9 NDS 3 + 0.0000754xl0" 12 NDS 4 


For gear 3: 

SPIN1 = -0.05786 + 0.6550xl0- 3 NDS + 0.02488xl0- 6 NDS 2 

+ 0.06789x10-9 NDS 3 _ 0.007325xl0‘ 12 NDS 4 

SPIN2 = -0.000167 + 0.008048x10-3 NDS - 0.000706xl0" 6 NDS 2 

+ 0.00193x10-9 nos 3 - 0.0001 554x1 O" 1 2 NDS 4 


For gear 4: 

SPIN1 = -0.03036 +0.4248x10-3 NDS + 0 . 0214x1 0~ 6 NDS 2 

+ 0.01863x10-9 NDS 3 - 0.001192xl0" 12 NOS 4 

SPIN2 = -0.0010137 + 0.01585x10-3 NDS - 0.005401 xlO" 6 NDS 2 

+ 0.001154x10-9 NDS 3 - 0. 00005956x1 O" 1 2 NDS 4 


For NDS < 10, HPLOSS = 0 


XII. Axle Efficiency 

EFAXL = 0.0785 + 0.505624CR - 0.098384CR 2 + 0. 00648856CR 3 
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rn ,„„/5252 HPDs\ 
CR ■ ,os \ NWHL j 

0 < CR < 5.983 

HPDS > O.lxlO" 4 

XIII. Wheel and Drive Shaft 


TQWHL = GRAXL(EFAXL)(TQDS) = 


5252 HPWHL 
NWHL 


TQDS = 525 m1c PDS ^ 


HPDS = HPOTR - HPLOSS 


HPWHL = = HPDS(EFAXL) 


XIV. Road Load 


HPRL 



[1.0 + ROLl(VEL)] + 0.002557 (<j) ( CD) ( AFNT)( VEL) 



+ WT( GRADE) 



XV. Vehicle Velocity 


dVEL _ 21 .937( FACC) 
dt “ WTEQ 


FACC 


TQACC 

RADIUS 


WTEQ = WT + 


322W3 

1 2( RADIUS) 2 


RADIUS = 


5280 

2ir REVMI 


_ 52 52j HPACC j 

IVHV.V. - MI.IUI 


HPACC = HPWHL - HPRL 
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1. Advanced Gas Turbine (AGT) Powertrain System Development for Automotive 
Applications. (Garrett No. 31-3725(3), Garrett Turbine Engine Co.; NASA 
Contract DEN3-167.) DOE/NASA/Ol 67-811 3, NASA CR-167901, 1981. 

2. Bruton, William M.; and Cwynar, David S.: Lewis Hybrid Computing System- 

Users Manual. NASA TM-79111, 1979. 
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TABLE I. - COMPARISON OF STEADY STATE HYBRID AND 
DIGITAL SIMULATION RESULTS 


(a) Throttle setting, 20° (b) Throttle setting, 40° (c) Throttle setting, 50° 


Engine 

vari- 

ables 

Digital 

simu- 

lation 

Hybrid 

simu- 

lation 

GAMMA 

50° 

50° 

IGV 

39.8 

39.9 

BETA 

50.22 

50.39 

GEAR 

4 

4 

NEP 

76.04 

75.92 

VEL 

70.08 

70.07 

WAIRL 

0.4131 

0.4130 

WFL 

12.53 

11.31 

T4 

2903 

2856 

T5.1 

2468 

2427 

T3.6L 

2401 

2374 

TQCAR 

2.620 

2.541 

TQACC 

12.106 

2.200 

TQDS 

60.42 

60.00 

TQOTR 

95.22 

90.00 

NOTR 

1897 

1896 

HPRL 

24.84 

29.37 

TQWHL 

158.2 

160.0 

TCTR 

1 .203 

1 .198 

TCKI 

373.2 

380.0 

1CSR 

0.7851 

0.7901 

NCTR 

2416 

2402 

EFAXL 

0.9590 

0.9589 


Engine 

vari- 

ables 

Digital 

simu- 

lation 

Hybrid 

simu- 

lation 

GAMMA 

40° 

40° 

IGV 

41 .67 

41 .62 

BETA 

74.82 

78.04 

GEAR 

4 

4 

NEP | 

67.45 

67.36 

VEL 

49.85 

49.84 

WAIRL 

0.3240 

0.3231 

WFL 

7.325 

6.985 

T4 

2769 

2792 

T5.1 

2437 

2457 

T3.6L 

2389 

2413 

TQCAR 

1 .649 

1 .610 

TQACC 

12.26 

11.00 

TQDS 

38.05 

30.00 

TQOTR 

61 .22 

54.00 

NOTR 

1349 

1348 

HPRL 

11.78 

12.93 

TQWHL 

99.16 

80.00 

TCTR 

1.379 

1 .325 

TCKI 

492.4 

530.0 

TCSR 

0.5573 

0.5603 

NCTR 

2420 

2406 

EFAXL 

0.9547 

0.9505 


Engine 

vari- 

ables 

Digital 

simu- 

lation 

Hybrid 

simu- 

lation 

GAMMA 

20° 

20° 

IGV 

45.00 

45.02 

BETA 

69.01 

68.18 

GEAR 

3 

3 

NEP 

63.00 

62.93 

VEL 

41 .30 

41 .32 

WAIRL 

0.2756 

0.2737 

WFL 

5.008 

4.050 

T4 

2688 

2716 

T5.1 

2417 

2442 

T3.6L 

2380 

2425 

TQCAR 

1 .144 

1 .148 

TQACC 

15.55 

15.40 

TQOS 

32.10 

35.10 

TQOTR 

35.06 

35.10 

NOTR 

1676 

1675 

HPRL 

7.627 

8.635 

TQWHL 

83.46 

92.00 

TCTR 

1 .024 

1 .017 

TCKI 

487.1 

485.0 

TCSR 

0.8752 

0.8797 

NCTR 

1915 

1904 

EFAXL 

0.9524 

0.9529 
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Figure 1. - Advanced gas turbine powertrain 
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WAIRL/1.0 VEL/100 NEP/110 GRTR/2.4 BETA/100 






WFL/50 WAIRL/1.0 VEL/100 NEP/110 GRTR/2.4 BETA/100 IGV/60 GAM/60 



VSTC vane response, deg GAMMA sfep command> deg 



(a) GAMMA step command. (b) Compressor inlet guide vane 

position response. 



o 10 20 30 40 50 60 70 80 " 0 10 20 30 40 50 60 70 80 


Time, sec 

(c) VSTC vane response. (d) Percent engine speed response. 

Figure 5. - Digital simulation dynamic results ; GAMMA step command from 0° to 2(P. 


Fuel flow response, Ib/hr V e hic ,e speed response, mph 



Response of driveshaft torque, ft-lb Response of torque at input of gearbox, ft-lb 



.38 

.36 

o 

S .34 

-O 

«? .32 


a .30 
»/> 



Time, sec 


(h) Engine torque response. 
Figure 6. - Continued. 


(g) Fuel flow response. 
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